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ABSTRACT 

The well-established correlations between the mass of massive black holes (BHs) in the nuclei of most stud- 
ied galaxies and various global properties of their hosting galaxy lend support to the idea that dwarf galaxies and 
globular clusters could also host a BH in their centers. Direct kinematic detection of BHs in dwarf spheroidal 
(dSph) galaxies are seriously hindered by the small number of stars inside the gravitational influence region of 
the BH. The aim of this Letter is to establish an upper dynamical limit on the mass of the putative BH in the 
Ursa Minor (UMi) dSph galaxy. We present direct N-body simulations of the tidal disruption of the dynamical 
fossil observed in UMi, with and without a massive BH. We find that the observed substructure is incompatible 
with the presence of a massive BH of (2-3) x 10^ Mq within the core of UMi. These limits are consistent 
with the extrapolation of the M^h - cr relation to the M^h < 10^ Mq regime. We also show that the BH may be 
off-center with respect to the center of symmetry of the whole galaxy. 

Subject headings: galaxies: individual (Ursa Minor dwarf spheroidal) — galaxies: kinematics and dynamics 
— galaxies: dwarf — stellar dynamics 



L INTRODUCTION 

Intermediate-mass black holes (IMBH; Mbh ^ 10^ - 
10'*Mq) accreting gas from their surroundings have been pos- 
tulated to explain the engines behind ultraluminous X-ray 
sources recently discovered in nearby galaxies (see Colbert & 
Miller 2005, for a review). IMBHs would fill the existing gap 
between stellar and supermassive black holes found in active 
galactic nuclei. According to the tight relation between Mbh 
and the central velocity dispersion a (Gebhardt et al. 2000; 
Ferrarese & Merrit 2000; Tremaine et al. 2002), or between 
Mbh and the mass of the bulge (Magorrian et al. 1998), a nat- 
ural place to look for these IMBHs are astrophysical systems 
less massive than normal galaxies, such as dense star clus- 
ters, globular clusters, and dwarf galaxies. If Mbh is cor- 
related with the total gravitational mass of their host galaxy 
(Ferrarese 2002; Baes et al. 2003), central BHs should be an 
essential element in dark matter dominated objects -such as 
dwarf spheroidal (dSph) galaxies. Estimates of the mass of 
BH in these systems would be of great importance in com- 
pleting the Mbh - c relation. 

Direct and indirect searches for IMBH at the centers of 
globular clusters and small galaxies have been attempted (e.g., 
Gerssen et al. 2002, 2003; Valluri et al. 2005; Maccarone et 
al. 2005; Ghosh et al. 2006; Maccarone et al. 2007; Noyola 
et al. 2008). For instance, Noyola et al. (2008) have reported 
a central density cusp and higher velocity dispersions in their 
central field of the globular cluster uj Centauri, which could 
be due to a central BH {Mbh ~ 4 x IQ'^Mq). However, a new 
analysis by Anderson & van der Marel (2009) and van der 
Marel & Anderson (2009) does not confirm the arguments 
given by Noyola et al. (2008), and provides an upper limit 
for the BH mass of L2 x IO'^Mq. 

There is sparse evidence that BHs could be present in at 
least some dSph. The Ursa Minor (UMi) dSph galaxy has 
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been also suspected to contain a BH of ^ 10^ Mq (Strobel 
& Lake 1994; Demers et al. 1995). Maccarone et al. (2005) 
discuss the possibility that the radio source found near the 
core of UMi is, in fact, a BH with a mass ^ Kf'MQ. 

In this Letter, we examine the dynamical effects of puta- 
tive IMBHs in the core of dSphs on the very integrity of cold, 
long-lived substructure as that observed on the northeast side 
of the major axis of UMi. Although UMi has long been sus- 
pected of experiencing ongoing tidal disruption, regions with 
enhanced volume density and cold kinematics cannot be the 
result of tidal interactions (Kleyna et al. 2003, hereafter K03; 
Read et al. 2006; Sanchez-Salcedo & Lora 2007). This sug- 
gests that the secondary peak in UMi is a long-lived struc- 
ture, surviving in phase-space because the underlying gravita- 
tional potential is close to harmonic (K03). Sanchez-Salcedo 
& Lora (2007) derived an upper limit on the mass and abun- 
dance of massive dark objects in the halo of UMi to avoid a 
quick destruction of the clump by the continuous gravitational 
encounters with these objects. In this work, we will assume 
that the dark halo is comprised of a smooth distribution of 
elementary particles and then study the disintegration of the 
clump, placing contraints on the mass of a possible central 
IMBH in UMi. 

2. INITIAL CONDITIONS AND W-BODY SIMULATIONS 
2. 1. Ursa Minor and its clump 

Ursa Minor, located at a galactocentric distance of Rgc = 
76 ±4 kpc (Carrera et al. 2002; Bellazzini et al. 2002), is one 
of the most dark matter-dominated dSphs in the Local Group, 
with a central mass-to-light ratio M /L > 100 Mq/L© (e.g., 
Wilkinson et al. 2004). The measured central velocity disper- 
sion is 17 ±4 km s"' (Munoz et al. 2005) and the core radius 
of the stellar component along the semimajor axis is ~ 300 pc 
(Palma et al. 2003). UMi reveals several morphological pecu- 
liarities: (1) The shape of the inner isodensity contours of the 
surface density of stars appears to have a large ellipticity of 
0.54, (2) the highest density of stars is not found at the cen- 
ter of symmetry of the outer isodensity contours but instead is 
offset southwest of center, (3) the secondary density peak on 
the northeast of the major axis is kinematically cold. 

The secondary density peak has a Icr radius of ~ 1 .6' (~ 35 
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Fig. 1. — Snapshots at t = 2, 8 and 13 Gyr for Riji = 50 pc and Rcok = 510 pc without self-gravity (top panel) and with self-gravity (middle panel). For 
comparison, the snapshots for Rcore = 300 pc and with self-gravity are also shown (bottom panel). 



pc at a distance of 76 kpc) when fitted with a Gaussian pro- 
file. The stellar distribution that forms this density excess is 
elongated not along the major axis of the isodensity contours 
of the elliptized King model, but along a line at an intermedi- 
ate angle between the major and minor axes of these contours 
(Palma et al. 2003). The bend in the isodensity contours in- 
dicates that such clump is gravitationaUy unbound. Interest- 
ingly, K03 found that the velocity of stars within a 6' (130 pc) 
radius aperture on the clump are best fitted by a two-Gaussian 
population, one representing the underlying 8.8 km s~' pop- 
ulation, and the other with a line-of-sight velocity dispersion 
of 0.5 km s"' . Although the value of the cold population is ill 
determined, we can be certain that the velocity dispersion is 
< 2.5 km s"' at a 95% confidence level. The mean velocity of 
the cold population is equal to the systemic velocity of UMi, 
implying that either the orbit is radial and the clump is now at 
apocenter or the orbit lies in the plane of the sky (K03). 



2.2. Initial conditions 



We consider the evolution of a clump inside a rigid halo of 
dark matter with a density law: 



p{r) = 



1/2' 



(1) 



where po is the central density and Rqok is the dark halo core 
radius. This profile was chosen in order to facilitate compari- 
son with K03. We explore different values for /?core- Follow- 
ing K03, once /?core is fixed, we rescale the central density to 
have a dark matter mass of 5 x 10'' Mq inside 600 pc, which is 
approximately the maximum extent of the stellar distribution 
observed in UMi. Therefore, in all our models, we have a total 
mass-to-Hght ratio of (M,o,/Lv) « 90Mq/Lq inside 600 pc, 
for a visual luminosity Ly = 5.4 x 10^ Lq (Palma et al. 2003). 

For a normal stellar population with M/Ly = 2Mq/L}^, the 
stellar mass within the core radius of the dwarf is ~ 4.5 x 
IO^Mq, whereas the dark matter mass is > 6.2 x lO^M©. 
Therefore, the contribution to the potential of the baryonic 
mass was ignored. 

The initial density profile of the clump follows a Plummer 
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mass distribution, 



p{r): 



M,.Ri 



(2) 



where = 2 x 10 Mq is the total mass of the clump and Rp 
is the Plummer radius. One should note that there is a simple 
relation for a Plummer model between Rp and the half-mass 
\3Rp. We use /?i/2 values between 25 and 



radius: R 



1/2 



50 pc to initialize our simulations. The clump's self-gravity 
is included to have a realistic and complete description of its 
internal dynamics because, as we will see later, the tidal radius 
of the clump may be larger than /?i/2 when large values of 
/?core are used. 

The clump is dropped at the apogalactocentric distance of 
200 pc from the UMi center with a certain tangential veloc- 
ity \'T, which defines the eccentricity e of the orbit. The orbit 
of the clump lies in the x-y plane, which is also the plane of 
the sky. For simulations with a central BH, it is clear that the 
evolution of the clump depends on its orbital eccentricity. For 
instance, we found that for a radial orbit, the BH dissolves 
the group of stars in its first passage through the galactic cen- 
ter. In order to provide an upper limit on the mass of the BH 
and since there is more phase space available for nearly cir- 
cular orbits than for radial ones, we will take, in most of our 
simulations, a rather circular orbit with e = 0.5, which corre- 
sponds to vj = 5.7 km s"' for Rcok =510 pc. However, to our 
surprise, the disintegration time was found to be insensitive to 
the initial eccentricity of the clump's orbit as long as e < 0.87. 
Altogether, there are essentially three model parameters that 
we have explored in our simulations: R^ok, Ri/i and Mbh- 

2.3. N-body simulations 

We developed an N-body code that links an individual 
timestep to each particle in the simulation. Only for the par- 
ticle that has the minimum associated time, the equations of 
motion are integrated (with a second order predictor-corrector 
method). This "multi-timestep" method reduces the typical 
CPU time of direct, particle-particle integrations (oc A^^, with 

the number of particles), and allows integrations of systems 
of '^lOOO particles to be carried out with relatively short CPU 
times. 

The cluster density is so tiny that the internal relaxation 
timescale is very long (~ 11 Gyr for an initial cluster with 
Ri/2 = 50 pc). Therefore, although our code is suitable to in- 
clude two- and three-body encounters, the cluster behaves as 
collisionless and internal processes such as evaporation do not 
contribute to the dissolution of the cluster. For the same rea- 
son, all the particles have the same mass, and the presence of 
binaries was ignored. 

All the simulations presented in this paper used 600 parti- 
cles, each one having a mass of 33M0. We chose a smoothing 
length of 0.7 pc, which is approximately 1/lOth the typical 
separation among cluster particles within 7? 1/2 at the begin- 
ning of the simulation. The convergence of the results was 
tested by comparing runs with different softening length and 
A^. The effect of adopting a dfferent smoothing radius be- 
tween 0.1 and 1 times the typical distance was found to be 
insignificant. We run the same simulation with different 
(A^ = 200,400,600 and 1800 particles), and convergence was 
found for A^> 400. 

In order to validate our simulations, we checked that, when 
the cluster evolves at isolation, the Plummer configuration is 
stationary and that the energy is conserved over 12 Gyr. To be 



certain that the interaction with the BH is resolved well, we 
compared the change in kinetic energy of the cluster, when 
colliding with a particle of mass 1O^M0, moving at 200 km 
s"', with the predictions in the impulse approximation (Bin- 
ney & Tremaine 2008), for different impact parameters, and 
found good agreement (differences less than 10%) between 
them. 

We were also able to reproduce K03 simulations of the 
evolution of an unbound clump (ignoring self-gravity) in the 
gravitational potential created by the mass distribution given 
by Eq. ([T]i- We confirm the K03 claim that even if the clump 
is initially very compact, with a Icr radius of 12 pc, cusped 
halos cannot explain the survival of the substructure for more 
than 1 Gyr, while the substructure can persist for ^ 12 Gyr in 
halos with core radii 2-3 times the clump's orbit. 

3. RESULTS 

3.1. Simulations without BH: the role of self- gravity 

We have first examined the evolution of a clump with 
^1/2 = 50 pc embedded in a galaxy with a large core of 5 10 pc 
(see Fig.[T]i. Note that when self-gravity is included, an initial 
Rij2 value of ^ 15 pc, as that used in K03, is no longer real- 
istic because the clump remains too compact over its lifetime 
to explain its present appearance. With self-gravity, strong 
substructure continues to persist for a Hubble time, while in 
the non-self-gravitating case the substructure is completely 
erased in 10 Gyr. The evolution of the same self-gravitating 
clump but in a galaxy with R^ok = 300 pc is also shown in Fig. 
[U As expected, when the size of the galaxy core is decreased, 
the clump does not survive so long, but we can have a galaxy 
core of 300 pc without causing the clump to disintegrate in 
13 Gyr. We confirm the K03 claim that cored halos allow the 
structures to remain uncorrupted for a Hubble time. More- 
over, we find that, when including self-gravity, a galaxy core 
of 1 .5 times the clump's orbit is enough to preserve the in- 
tegrity of the clump. 

3.2. Simulations with BH 

If UMi hosts an IMBH, the density substructure is erased 
not only by the orbital phase mixing and the tidal field of the 
galaxy but also by the gravitational interaction with the hypo- 
thetical BH. Since the destructive effects depend on the mass 
of the BH, we can establish an upper limit on the mass of 
the BH in UMi by imposing that the BH must preserve the 
longevity of the clump for more than 10 Gyr. We proceeded 
to add a massive particle of 3 x 10"* Mq at the center of the 
galaxy potential, emulating a BH. Figure |2] shows the evolu- 
tion of a self-gravitating clump with an initial size of 7J 1/2 = 50 
pc in a galaxy with /?core = 510 pc. Under the influence of the 
gravitational pull exerted by the clump, the BH, initially at 
rest, is displaced from the center of UMi. The azimuthal or- 
bital lag angle of the BH in the {x^y) plane is w 7r/2. As a 
consequence, the clump feels a gravitational drag and loses 
angular momentum which is transferred to the BH. Due to 
the angular momentum loss, the clump starts to spiral into 
the center. At 2 Gyr, we can see that, in fact, the clump has 
smaller orbital radius than it would have in the absence of 
the BH (see Fig. |2]l. The BH, on the other hand, spirals out- 
wards to larger radii until it reaches its maximum galactocen- 
tric radius. At f = 4 Gyr, the clump reaches the galactic center 
and this "exchange" of orbits starts all over again. Therefore, 
when looking at the BH of UMi, one should bear in mind 
that the BH does not necessarily settle into the center of the 
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Fig. 2. — Snapshots at 2, 8 and 13 Gyr for a simulation with R|^2 - 
mass Mbh = 3x10'' Mq (bottom panel). 



50 pc and Score = 5 10 pc (top panel), and the same simulation but with a BH (triangle) with 



galaxy. It would be interesting to perform N-body simula- 
tions including the stellar background to elucidate if the two 
observed off-centered regions with the highest stellar density 
(e.g., Palma et al. 2003) are a consequence of the dynamical 
response to the BH. The minimum distance between the BH 
and the center of mass of the clump is r-.^ 100 pc in this model. 
Since the orbits of the clump and the BH never cross, the tidal 
disruption of the group of stars can be described as a secular 
process of stellar diffusion in phase-space. This slow relax- 
ation process forms a stellar debris of stripped stars that move 
on a galactic orbit very similar to that of the clump itself. 

In Fig. |2]we see that, under the influence of the BH, the 
clump is completely dissolved at 8 Gyr. In order to quan- 
tify the clump's evolution, we calculated a map of the sur- 
face density of particles in the ix,y) plane at any given time t. 
We sample this two-dimensional map searching for the parcel 
(of 20 X 20 pc size) that contains the highest number of stars 
A/'max(0- This region is centered at the remnant of the clump. 
The number of stars A/lnax as a function of time is shown in 
Fig. [3] As expected, A/"max decreases as the simulation evolves 
due to the spatial dilution of the clump caused by tidal heat- 
ing. Once A/^ax drops a factor 2, which occurs in r^ 6 Gyr, the 
tidal disruption process is accelerated and Afmux dramatically 
declines at ~ 8 Gyr (around the orbit 55), which can be taken 
as the disruption time, denoted by f^. The disruption of the 
clump was confirmed by visual inspection of the simulations. 

Further simulations show that even if the eccentricity of the 
orbit is similar to that of the isodensity contours of the sur- 
face density of UMi background stars, the disruption time is 
essentially the same. 

We carried out simulations with /?i/2 = 50 pc, Rcore = 510 
pc and e = 0.5, but with different Mbh- We found that the 
substructure disintegrates in 11 Gyr for a 10"* Mq BH, and in 
2 Gyr for Mbh =10^ Mq. We also explored different com- 



binations for Rrnre and R 



1/2- 



A reduction of the core of the 



galaxy leads to a more stringent upper limit on the mass of 
the BH. However, one can increase the longevity of the struc- 
ture if a smaller value for 7? 1/2 is adopted. For instance, we 
also find that frf ~ 8 Gyr for /^core = 200 pc, R1/2 = 25 pc and 
Mbh = 3x IO^Mq. 

In the absence of any knowledge about the initial dynamical 
state of the BH, the most natural assumption is that the center 
of the galaxy was its presumed birth site^. For completeness, 
we have also carried out simulations for an off-centered BH. 
A compilation of the models is given in Table[T] When the BH 
is on radial orbit, ~ 4-9 Gyr, for Mbh = 3 x W^Mq. Only 
in very special circunstances -when the BH is not on a radial 
orbit, the orbits are coplanar and the azimuthal lag angle is 
TT-, the disruption of the cluster is less efficient because the 
separation between the clump and the BH is larger on average 
over the simulation. 

We conclude that, if UMi's dumpiness is a primordial arti- 
fact, then the survival of the secondary peak imposes an upper 
Umit on the mass of the putative BH of Mbh = (2-3) x 10'* 
Mq, if the BH originally lurked at the center The extrapola- 
tion of the Mbh - c relation for elliptical galaxies (Gultekin et 
al. 2009) predicts a 1.0 i^ *^ x 10"* Mq BH for UMi. There- 
fore, our constraint is still consistent with both the extrapo- 
lated value and that inferred by Maccarone et al. (2005). 

4. DISCUSSION AND CONCLUSIONS 

Our self-gravitating simulations confirm the claim of K03 
that the dark halo of UMi must have a core radius of ^ 300 pc, 
comparable to the core radius of the underlying stellar popu- 
lation, in order to preserve the integrity of this clump. While 

^ The BH candidate in UMi reported by Maccarone et al. (2005) is placed 
at about 7 pc from the center. 
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TABLE 1 

Relevant parameters and destruction times. 



Run 




«l/2 


BH position* 


BH velocity 


Mbh 


td 


# 


pc 


pc 


at r = [pc] 


at f = [km s"'] 


Me 


Gyr 





510 


50 


— 


— 





> 14 


1 


510 


50 








5 X 10^ 


> 14 


2 


510 


50 








1 X 10" 


11 


3 


510 


50 








3 X 10" 


8.3 


4 


510 


50 








1 X 10' 


1.7 


5 


200 


25 








3 X 10" 


7.8 


6 


510 


50 


(0,-50,0) 





3 X 10" 


5.5 


7 


510 


50 


(0,0,-50) 





3 X 10" 


8.5 


8 


510 


50 


(0,-50,0) 





5 X 10" 


3.5 


9 


510 


50 


(0,0,-50) 





5 X 10" 


5.0 


10 


510 


50 


(0,0,-50) 





1 X 10' 


1.8 



20 



* The clump is initially at the position (200,0,0). 

orbits 

40 60 80 



1 00 



10 - 




Fig . 3 . — Smoothed curves of Afmm- as a function of time for runs # (black 
line), # 3 (grey Hne), # 5 (red line) and # 6 (blue line). 

the firm dynamical detection of an IMBH in any dSph galaxy 
is challenging because it requires observations of the veloc- 
ity dispersion of stars deep into the core, we have demon- 
strated that the very integrity of kinematically cold substruc- 
ture in dSph galaxies may impose useful limits not only on 
the core of the galaxy but also on the mass of putative BHs. 
In the case of UMi, the maximum mass of a BH initially 
seated at the centre of the potential or initially on radial or- 
bit, is (2-3) X 10"* Mq. When searching for direct detection 



of the possible IMBH in UMi (e.g., Maccarone et al. 2005), 
one should keep in mind that the BH may be offset from the 
galactic center because of the gravitational pull exerted by the 
clump. 

As pointed out by the referee, the contour map of the sur- 
face brightness of the nucleus of M3 1 also has a bright off- 
center source (Fig. 2 of Lauer et al. 1993), which was inter- 
preted not as a separate stellar system but as the apoapsis re- 
gion of an eccentric stellar disk orbiting a central massive BH 
(Tremaine 1995; Lauer et al. 1996). The timing problem of 
the clump may be circumvented if one assumes that the sec- 
ondary peak of UMi is the part of a ring close to apoapsis. 
Since the stars forming the secondary peak in UMi have a 
mean velocity equal to the systemic velocity of UMi, the ring 
should lie close to the plane of the sky and should be very 
eccentric. A serious difficulty with this scenario is that the 
apsides of the orbits need to be extremely aligned. Our sim- 
ulations have shown that this alignment is not possible if the 
ring is the tidal debris of a stellar cluster. We do not know 
a satisfactory mechanism to explain the formation of a very 
eccentric stellar ring at scales of a galaxy with the required 
apsidal alignment. 
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